ON THE FRACTIONAL BLACK-SCHOLES MARKET 
WITH TRANSACTION COSTS 

EHSAN AZMOODEH 

Abstract. We consider fractional Black-Scholes market with propor- 
tional transaction costs. When transaction costs are present, one trades 
periodically i.e. we have the discrete trading with equidistance n~^ be- 
tween trading times. We derive a non trivial hedging error for a class of 
11^ I European options with convex payoff in the case when the transaction 

costs coefBcients decrease as n~^^~^' . We study the expected hedging 
jO' error and asymptotic behavior of the hedge as H —>■ h. 
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Absence of arbitrage is a weak form of financial equilibrium and it plays a 
fundamental role in mathematical finance. In fractional Black-Scholes mar- 
ket, stock price is modelled by geometric fractional Brownian motion. This 
model admits arbitrage opportunities without transaction costs with con- 
tinuous trading. Guasoni [10] showed that, if one considers such financial 
markets with proportional transaction costs, then arbitrage opportunities 
disappear. In [IjJ, the authors studied a more general formulation of such 
financial markets with proportional transaction costs and they also find a 



t*"^ I superreplication price for a class of European options. 



Leland [15] suggests a way to include proportional transaction costs in clas- 
sical Black-Scholes market. Namely, one trades in discrete and equidistant 
trading times i.e. periodically. Over each trading subinterval, the trader 
cH I follows the delta hedging strategy computed at the left point of the trading 

subinterval with a modified volatility. The modified volatility depends on 
the original volatility and the number of trading intervals (see |14)). Leland 
remarked that the price of the modified strategy approximately hedges the 
option payoff at terminal date as the length of trading intervals tends to zero. 

The Leland's approach is considered useful for practitioners, although the 
modified strategy does not provide an exact hedge in the case when the 
level of transaction costs is constant. First Lott [IS] and later Kabanov and 
Safarian |14] showed that this is true when transaction cost coefficients kn 
decrease as 

kn = kon~°', ae(0, -], 
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where - is the length of the trading intervals. Constant transaction costs 
correspond to the case a = 0. Kabanov and Safarian [H] computed the 
hedging error explicitly for European call option and in general it is not 
equal to zero. They showed that the limiting hedging error is positive with 
probability one. This means that the option is always underpriced in the 
limit. 

The motivation of this study comes from the recent work by Azmoodeh et. 
al. [2]. There we studied a hedging problem for European options with con- 
vex payoff in fractional Black-Scholes market with Hurst parameter H > 2- 
We assume that the market is frictionless and continuous trading is possi- 
ble. Then any European option with convex payoff can be hedged perfectly. 
Moreover, hedging strategy and hedging cost (see [3]) are given explicitly. 
Simply speaking we showed that the classical chain rule holds for convex 
functionals of geometric fractional Brownian motion. Moreover, the wealth 
process, which is stochastic integral of replication portfolio with respect to 
stock price process S, is the limit of Riemann-Stieltjes sums almost surely. 
This makes our model more interesting from financial point of view (see [B] 
and ^). 

European call option with strike price K serves as a motivating example for 
us. In this case the stop-loss-start-gain strategy 

ut = l{5t>A'} 

is self-financing replication strategy in our model. Note that this replicat- 
ing strategy is not self-financing in the classical Black-Scholes market with 
standard Brownian motion. This strategy is of unbounded variation, and 
hence it is not practical for our model with proportional transaction costs 
and continuous trading. One possibility is to trade periodically, and the 
level of transaction costs is the function 

of the number of trading intervals. This is similar to Leland [T5] in the case 
of classical Black-Scholes model with transaction costs. 

The paper is organized as follows. Section 2 includes some auxiliary results 
to handle the transaction costs term in the limit. Also the description of 
the model in precise way is given. Section 3 contains the main result. The 
paper ends with more discussion on specific example European call option 
and conclusion in the section 4. 



2. Preliminaries 

Throughout the paper (fi, T, P) stands for a complete probability space. 
Assume B^ = {Bf^)i^]n^T] be a standard fractional Brownian motion with 
Hurst parameter H G (0, 1) and St = ^o expjS/^} be geometric fractional 
Brownian motion, 5*0 € M+. 
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2.1. Auxiliary results on convex functions. We recall some results on 
convex functions. First, recall that every convex function / : M — )• M has a 
left-derivative f'_ and a right-derivative /^. 

The next theorem gives information about the left-derivative f'_ and right- 
derivative f\. 

Theorem 2.1. [19j The functions f_ and f_^ are increasing, respectively left 
and right- continuous and the set {x : /_(x) ^ f_^_{x)} is at most countable. 

Moreover, the second derivative of a convex function / exists as a distribu- 
tion, and first derivative can be represented in terms of the second derivative. 

Theorem 2.2. [19J The second derivative f of convex function f exists in 
the sense of distributions, and it is a positive Radon measure; conversely, for 
any Radon measure jjl on M, there is a convex function f such that / = /i 
and for any interval I and x G int(I) we have the equality 

(2.1) f'_{x) = ^ sgn(x - a)fj,{da) + a/, 

where aj is a constant and sgn x = 1 if x > and —1 if x <0. 

Remark 2.2.1. If the supp(fj,) is compact, then one can globally state that 

(2.2) fUx) = ^ j sgnix- a) fi{da) 
up to a constant term. 

2.2. Linear approximation of convex functions. Let / : M — t- M be a 
convex function. For each interval [a,b], let vr = {a = ag < oi < ... < a„ = 
6} be a partition of the interval and 

||7r|| := max (oj — Oj-i). 

l<i<n 

A piecewise linear function through points (aj,/(aj)) is called a convex lin- 
ear approximation of convex function / on the interval [a, b] based on the 
partition tt. 

Theorem 2.3. Let [a, b] be a closed interval and {iTm} be a sequence of 
partitions of the interval [a, b] such that 

IKmll —^0 as m, — ;■ oo 

where -Km = {ai,a2, ...,a„(m)}. 

Let Pm be a convex linear approximation of convex function f on the interval 
[a,b] based on partitions Tim- Then we have 
(i) On the interval (a, b) as m ^ oo 

Pm -^ f o,nd {Pm)- -^ /- pointwise. 

(ii) For any bounded continuous function g we have 



gdfim — ?■ / gdfj, as m — )• oo, 

a,b] ■'[i,^] 
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where jim stands for Radon measure corresponding to the second de- 
rivative of Pm- 



2.3. Approximation of local time of fractional Brownian motion. 

The occupation measure related to fractional Brownian motion is defined by 



V(^ xU) = X{tGl:Bl' GU}=J^ l^B^eU} 



dt 



where / and U are Borel sets on time interval [0, T] and the real line re- 
spectively and A stands for Lebesgue measure. It is well-known that the 
occupation measure has a jointly continuous density {local time) which is 
denoted by l^{x, t) := l^{x, [0, t]) and is Holder continuous in t of any order 
a < 1 — H and in x of any order /3 < ^^ir (^°^ ^ survey article on the subject 
see Geman and Horowitz [9j). 

Let B^ = {-B^(i)}tg[o,T] be polygonal approximation of size A of B^ i.e. 
B^ is the polygonal lines which connect points {{iA, B^ (iA))} for suitable 
running index i. Set 

C^iB", [0, T]) = {te [0, T] : B^{t) = a and t ^ iA for each index i} 

and 

7vx(i?^,[o,r]) = #Q(i?^,[o,r]), 

i.e. the number of level a crossing of B^ over interval [0,T]. Then we have 
the following approximation for the local time l^{a,t). 

Theorem 2.4. Assume B be a standard fractional Brownian motion with 
H G (0, 1) and N^{B^ , [0, T]) be the number of level a crossing of size A— 
polygonal approximation of B^ . Then 



^A'-"NliB",[0,T])^l^ia,[0, 



T]) in L^ as A^ 0. 



Proof. See [T], Theorem 5. 



2.4. The model. Let / : M — t- M be a convex function with positive Radon 
measure /i corresponding to its second derivative. We showed in [2\ that 
in a geometric fractional Brownian motion market model with H > ^^ sll 
European options with the terminal payoff fiSx) can be hedged perfectly 
and the hedging strategy and hedging cost are f^{St) and f{So) respectively. 
More precisely we have that 



(2.3) f(^ST) = fiSo)+ I fUSt)dSt. 

Jo 

The stochastic integral in the right-hand side is understood as a limit of 

Riemann-Stieltjes sums almost surely, i.e. 

(2.4) J2 fUs^_,){s^ - s^_^) ^ f fUst)dSt, tr = f • 
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Assume that the terminal trading time T = 1. For each n we divide the 
trading interval [0, 1] to n subintervals [i"_i,if] where 

t1 = -=iAn, i = 0,l,...,n. 
n 

The two-assets market model consists of : 

(i) Riskless asset (bond), Bt = I; t G [0, 1] which corresponds to zero 
interest rate. 

(ii) Risky asset (stock) whose price is modeled by geometric fractional 
Brownian motion 

St = Soe^"; t G [0, 1], and H > ^. 

Now consider the discretized version of hedging strategy, i.e. 

n 
i=l 

In the presence of proportional transaction costs, the value of this portfolio 
at terminal date is 

(2.5) Vi(r) = f{So) + f e^dSt - kJ2St- j£(%) - fUSt- ,)|. 

Here it is assumed that the transaction costs are "two-sided" i.e. buying and 
selling are equally charged and transaction costs coefficient A; is a function 
of the length of trading intervals, i.e. 

k = kn = kon^'^^'^K ko>0. 

Remark 2.4.1. In general the hedging strategy f_{St) is not of bounded 
variation. The strategy 

f^{St) = l{St>K} 
is of unbounded variation, for f{x) = {x — K)'^ corresponding to European 
call option with strike price K . 

Remark 2.4.2. Note that it is assumed that there is no transaction costs 
at time t = when the trader enters to the market. 

3. Main Result 

After the preliminaries and the precise description of the market model, we 
can state the main result of the paper. 

Theorem 3.1. Assume the level of transaction costs k = kn = kQn~'^^~^' , 
where ko > 0. Then 

p- hm yi(r) = /(^i)-J, 

n— >oo 

where 

J = J(A;o) := J-fco I I Stl" {In a, dt)n{da), 
V vr Jji^Jq 
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and the inner integral in the right hand side is understood as limit of Riemann- 
Stieltjes sums a.s. 

When the number of portfoho revision increases fast enough or equally say- 
ing, the transaction costs coefficients kn decrease faster, we have the perfect 
replication in the limit, i.e. we have the following result. 

Corollary 3.1. Let the level of transaction costs k = kn = kon'" where 
a > 1 — H , kQ > 0. Then 

p- hm yi(r) = f{Si). 

Proof: We can assume that the support of fx is compact otherwise one can 
consider auxiliary convex functions 



fn{x) 
(see [2] for details). Since 



/;(0)x + /(0), ifx<0, 

f{x), ifO<x<n, 

f-{n){x — n) + f{n), if x > n. 



f{S,) = f{So)+ / f^{St)dSt, 



^i(^") = f{So) + E/-('^*?-i)('^*r - Sq_,) - kn^StnJf'_{Stn) - f'_{Se^_^^ 



i=l 



i=l 



Therefore, we have the indentity 



where 



v^[e^)-f{Si) = i\-koil 

II = Y.f'.{Svi i)(% - % J - / f'-{St)dSu 

n 

II = Ai-"^St^^JfUSt-) - fUStfJl 



i=l 



Note that Z^ — )• almost surely by ()2.4p . It remains to study the behavior 
of the second term I^. The representation ()2.2p relates the left derivative of 
convex function in term /^ to Radon measure of its second derivative. We 
divide the proof in three steps, depending on the supp /j,. 

Step 1: supp /i = {a}. 

We can assume that fi{a) = 1 and f_{x) = l{x>a}- This follows from 

representation (j2.2p . 
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For any m > n, 



m -1 

^m 2^5't"lJl{S'tm>a} - ll^tm >a}\ - / Sfl (In a, dt) < 
~{ i J-1 Jo 

m n 

^^ml^\^Strp \l{s,n.>a}-'^{S,m >a}\-^Str^ ^ |l{5jm>a}-l{Stm >a} 

I ■' .? .7-1 .? J-1 



JG^W 



+ 






i=l 






^ri,m + tin,m ~r C/ji, 



where for each i = 1, 2, .., n, /(i) = {j : t^ G (^"-1; *?]}■ 
Obviously, hni^^oo C'n = 0, since /^(Ina, •) is increasing in t and 



;^ 



Hi 



r (x, (s, t]) = /^ (x, t) - r [x, s), s<t. 



\Bn,m\ < 2j%_i ^m'^ zJ I l{5tm>a}-l{Stm >a} I "^^(In a(t"_i, t"]) 

7 j — 1 



i=l 



Therefore for each fixed n by Theorem 12.41 as m — )■ oo, 



P- hm 

m— >C!0 



^m ^ V] |l{5tm>a} - l{5tm >a}l " ^^(In a, (t"_i, t"]) 

7 7 — 1 



P- Hm Ai-^iVi-(i?^,(t" ,t"])-/^(lna,(C_i,tr]) 



Hence -Bn,m converges in probabihty to zero as m tends to infinity. 



\An,m\ <22^m^ z2 '^-i ~ 'S't™_i 1 1 l{Stm >a} " 1 



j=l 



je^W 



>a} -"-{Stm >a}l 



< yZ sup |5i-_^ - Su\Al^ ^ yZ |l{5im>a} - l{5im > 



^ineCtj'.i,*^) 



a} I 



i6/{i) 



A ^ sup \St^ ^ - 5„|/^(lna, (^1,*^]) = ^n- 

Fix e > 0. Since St is uniformly continuous on compact intervals [tf_;^,tf], 
there exists no G N such that for all n > no we have 

sup \St^_ — S'ul < e, 1 < i < n. 

«e(tr-i-tr) 
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Therefore, for all n > riQ, 



An < £ V /^(Ina, {tti,t:]) = e/^(lnG, [0, 1]). 



i=l 



So An converges to zero almost surely as n tends to infinity. 

Step 2: snpp fj, = {ai,a2, ■■■,ai}. 

Before to show the convergence in this case we need the following simple 
lemma. 

Lemma 3.1. For x,y gM. and positive numbers ai,a2, ■■■, ai we have 



^i'^{y>a,} - l{x>a,})ail = ^\'^{y>a,} - l{x>a,}|ai- 



Proof: First note that l{y>a.} — l{x>a } = ±1,0. Assume y > x, so this 
implies that all coefficients Isy^a } ~ ^{x>a } ^or a^-'s are equal to 1 or for 
1 < J < /. Hence 



I I 

I 



The case x > y is similar. 
Next we work with /^. 



'?!• 



i=l 
n ^ I 

I ^ z z z — i z — 1 

i=l j=l 

n -. I 

n ^ I 



i=l j=l 

An + -Dn 
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By lemma [STTj we see that 

^ n I 



i=l j=l 












Stf_ 


1 ^{St^>aj} l{Stn_paj} 




-^lyif^i^j) [ Stl"" (In a„dt) = 

^~i Jo 


5 li' «"<'"" 


dt)fi(da 


a similar argument for the term Bn we 


conclude that 




P- lim /2 = 

n->oo 


L 


Jo 


(In a, dt)/i(da). 





Step 3: general case. 

Let compact interval [a, b] contain the supp /x and let P^ be the convex 
linear approximation of convex function / on the interval [a, b] based on 
equidistant partition of interval [a,b] i.e. polygonal with vertices {(a + 
iAm{b — a), /(o + iAm{b — a)))}[^g. Define Pm's the same as / outside the 
interval [a,b]. Note that outside of the interval [a,b] the convex function / 
is linear. Then for m, n G IN we have 



\T'^ — tl < \T'^ — T l-t-lT" —T\-i-\T — Tl 



where 



In,m = ^n 2^ 'S't^-i I (-Pm)_ (5*4^) " {Pm)-{St^_^)\, 
i=\ 

Im= / Stl^ {In a, dt)^j.m{da). 
JnJo 
By the elementary inequality ||o| — |6|| < |a — 6|, we have that 

n 

j=l 
Now for fixed n, by Theorem 12 . 31 the right-hand side converges to zero almost 
surely as m tends to infinity. By simple calculations 

-1 /.I 



/ Stl"{lna,dt)= Stl{BH=ina}dt 



1 

Inai 




It follows that 



= al^(lno, [0, 1]). 
\Im - J\ < I / al" {In a, [0, 1]) fimida) - / al" {In a, [0,1]) fi{da) 

JR Jr 
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By Theorem 12.31 the right-hand side converges to zero almost surely as m 
tends to infinity, since support of /i is compact. 



Finally, for fixed m by Step 2, we have that 



P- lim I„jn = Im- 



Remark 3.1.1. Clearly, the hedging error 



J(fco) 



iH, 



ko / Sfl (In a, dt)fi{da) 
TT Jr Jo 

ko / a/^(lna, [0,l])/i(da) 
TT Jr 



is positive a.s. and strictly positive on a set of positive probability. So with 
proportional transaction costs, the discretized replication strategy asymptot- 
ically subordinates rather than replicate the value of convex European option 
f{Si) and the option is always subhedged in the limit. 

Remark 3.1.2. The limiting hedging error J = J(A;o) is small for small 
values of fixed proportional transaction costs coefficient /cq ■ 

4. Discussion and conclusion 

4.1. The case European call option. Consider European call option with 
corresponding convex function f{x) = {x — K)~^ . The approximating wealth 
process can not hedge perfectly the option payoff {Si — K)^ and limitting 
hedging error takes the form 

J = i^/^ (In is:, [0,1]). 

So, it is interesting to examine the expected hedging error 

IE(J) = i^]E(;^(lni^, [0,1])) 

^ /\-^exp{--t-2^1n2K}dt. 
Jo 2 



27r Jo 

The graph of the expected hedging error as a function with respect to vari- 
ables strike price K and Hurst parameter H is plotted in below which points 
out that the strike price K = 1 \s a, critical point. The expected hedging 
error tends to zero as strike price becomes bigger and bigger. 




Figure 1 . Expected hedging error 
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4.2. Asymptotic behavior with respect Hurst parameter H. Assume 
B^ be a standard fractional Brownian motion with Hurst parameter H G 
(0, 1). It is straightforward to check that when the Hurst parameter H tends 
to Ho E (0, 1), the finite-dimensional distributions of fractional Brownian 
motions {B } tend to finite-dimensional distributions of B °. It follows 
because of convergence of covariance functions and by 

JE\Bf - Sf |2 = \t- s|2^ 

and Billingsley criterion (see [5] ) we conclude that the family of laws of the 
fractional Brownian motions {B^} converge in law in the space C[0,r] of 
continuous function on the interval [0, T] to that of B^°. 

Moreover, the following theorem shows that the same holds for the family 
of local times {l^ , H G (0, 1)} of the fractional Brownian motions {B^}. 

Theorem 4.1. [13] The family {l^ ,H G (0, 1)} of local times of the frac- 
tional Brownian motions {B ,H G (0,1)} converges in law in the space 
C{[-D,D] X [0,r]), for any D,T > 0, to the local time l^° of fractional 
Brownian motion B^", when H tends to Hq. 

Now consider the family of fractional Brownian motions {B , i? G (2 , 1)}- 
Denote Brownian motion B and its local time l{x,t) as the limit in law when 
H tends to 2 of corresponding fractional Brownian motions {B } and its 
local times. Let S = S()exp{B} be geometric Brownian motion and f{S) 
be the payoff of European call option in Black-Scholes market model. Then 
we have that 

Corollary 4.1. Assume the fixed proportional transaction costs coefficient 
ko = \/\- Then 

lim ( lim T/(r)) = f(S)-Kli^uK,-), 

which limits take place in law in the space C[0,T]. 

Corollary 4.2. Assume the fixed proportional transaction costs coefficient 
ko = \\/\- Then 

lim ( lim l^(r)) = f(S) - Kl{lnK, •), 

which limits again take place in law in the space C[0,T]. 

4.3. Conclusion. It is known that geometric fractional Brownian motion 
with Hurst parameter H G (2;!) has zero quadratic variation process. 
Cheridito [7] uses this fact to show pricing model based on it, admits arbi- 
trage. Moreover in our frictionless model in the case of European call option 
the formula ()2.3p implies that 

{St - Ky = (So -K)++ [ l{s,>K}dSt. 

Jo 

This indicates that out-of-money options are worthless. Moreover the hedg- 
ing strategy ut = ^{St>K} is an arbitrage opportunity. On the other hand, 
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consider a frictionless pricing model with continuous price process 

Xt = So exp{Sf + sWt} e>0, 

where TV is a standard Brownian motion independent of B . This process 
has non-zero quadratic variation and fullfihs all conditions of a recent result 
by Bender et.al. [3]. Their result asserts that, this pricing model does not 
admit arbitrage opportunities with reasonable trading strategies. This in- 
dicates that the existence of non-zero quadratic variation is important for 
option pricing based on no arbitrage. 

Let / be a convex function with positive Radon measure fi as its second 
derivative. Consider continuous semimartingale Xt = X^e^* with local time 
Ix and Xq G M-|-. By Ito-Tanaka formula (see |19) . page 223) we have 

f{Xi) = f{Xo) + / f'_{Xt)dXt + \ f lx{a, [0, l])fi{da) 
Jo ^ JR 

= f{Xo) + / fUXt)dXt + \ I alw{\na, [0, l])^x{da). 

Jo ^ Jr 

Hence our fractional Black-Scholes model with asymptotic proportional trans- 
action costs has the same effect as the model which the stock price is mod- 
elled by semimartingale X. In this sense there is a connection between 
transaction costs and quadratic variation. 
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